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Abstract. In this paper we address the question whether 
age and metallicity effects can be disentangled with the 
aid of the broad-band colours and spectral indices from 
absorption feature strengths, so that the age of ellipti- 
cal galaxies can be inferred. The observational data un- 
der examination are the indices and [MgFe], and the 
velocity dispersion X for the sample of galaxies of Gonza- 
les (1993), supplemented by the ultra-violet data, i.e. the 
colour (1550- V), of Burstein et al. (1988). The analysis 
is performed with the aid of chemo-spectro-photometric 
models of elliptical galaxies with infall of primordial gas 
(aimed at simulating the collapse phase of galaxy forma- 
tion) and the occurrence of galactic winds. The galaxy 
models are from Tantalo et al. (1995). The study con- 
sists of four parts. In the first one, the aims are outlined 
and the key data are presented. In the second part, we 
summarize the main properties of the infall models that 
are relevant to our purposes. In the third part we present 
the detailed calculations of the spectral indices for single 
stellar populations and model galaxies. To this aim, we 
use the analytical relations of Worthey et al. (1994) who 
give index strengths as a function of stellar parameters. 
In the last part, we examine the age- metallicity problem. 
In contrast with previous interpretations of the Hp and 
[MgFe] data as a sort of age sequence (Gonzales 1993), 
we find that the situation is more complicate when the 
space of the four variables H^, [MgFe], (1550- V), and X 
is examined. Galaxies in the and [MgFe] plane do not 
follow a pure sequence either of age or metallicity. The 
observed (1550-V) colours are not compatible with young 
ages. Basically, all the galaxies in the sample are old ob- 
jects (say as old as 13-^15 Gyr) but have suffered from 
different histories of star formation. Specifically, it seems 
that some galaxies have exhausted the star forming activ- 
ity at very early epochs with no significant later episodes. 
Others have continued to form stars for long periods of 
time. This is perhaps sustained by the analysis of the gra- 
dients in the EL3 and [MgFe] indices across the galaxies. 
There are galaxies with no age difference among the vari- 



ous regions. There are other galaxies in which large gradi- 
ents in the mean age of the star forming activity between 
the central and the peripheral regions seem to exist. The 
nucleus turns out to be younger and more metal-rich than 
the outer regions. Finally, there are galaxies in which the 
nucleus is older but less metal-rich than the external re- 
gions. All this perhaps hints not only different histories 
of star formation but also different mechanisms of galaxy 
formation difficult to pin down at the present time. From 
the analysis of the H^, [MgFe], (1550-V), and X space, 
and of the age and metallicity gradients in single galax- 
ies, the suggestion is advanced that the overall duration of 
the star forming activity is inversely proportional to the 
velocity dispersion X (and perhaps galactic mass). 
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1. Introduction 

These days, it is understood that age and metallicity gov- 
ern the properties of the stellar populations in galaxies of 
different morphological type. Unfortunately, even in the 
simplest case of early-type galaxies (ellipticals), age effects 
often mask metallicity effects so that separating the two is 
a cumbersome affair (cf. Renzini 1986; Buzzoni et al. 1992, 
1993; Worthey 1994 for more details). The age-metallicity 
dilemma originates from the fact that increasing either 
the metallicity or the age makes the integrated spectral 
energy distribution (ISED) of a single stellar population 
(SSP) redder. 

The subject of galaxy ages has been tackled with differ- 
ent techniques of populations synthesis (evolutionary and 
optimizing, broad-band colours and narrow-band indices) 
with rather contrasting results. 

From evolutionary synthesis models and broad-band 
colours the indication arises that elliptical galaxies are 
made of old stellar populations moderately metal rich and 
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age of about 15 Gyr. Supplementary episodes of star for- 
mation at more recent epochs seem to be excluded. 

In contrast, the optimized synthesis models (cf. O' 
Connell 1986 and references) suggest that large age 
spreads are possible with a substantial contribution from 
young stars to the integrated light (for the prototype 
galaxy M32 as young as 5 Gyr). 

Recently, the study of spectral indices reveals a pow- 
erful tool to discriminate age from mctallicity effects. Ac- 
cording to Buzzoni et al. (1992, 1993) the use of calibrated 
Mg2 and Fig indices strengthens the conclusion reached 
with the broad-band colours. On the contrary, according 
to Gonzales (1993 and references) the indices [MgFe] and 
Hp hint that a substantial age spread ought to exist. 

The problem is further complicated by the poor knowl- 
edge of the star formation history of elliptical galaxies. 
Indeed the view that all elliptical galaxies are primordial 
old objects (cf. Bower et al. 1992a, b) coexists with the 
view that every elliptical galaxy is the product of merger 
events involving smaller galaxies (cf. Schweizer & Seitzer 
1992, Alfensleben & Gerhard 1994, Chariot & Silk 1994). 

In this work we would like to try a more articulated 
analysis aimed at separating age from metallicity effects. 
The study is based on the recent models for elliptical 
galaxies by Bressan et al. (1994), Tantalo (1994), and Tan- 
talo et al. (1995), however implemented with the calibra- 
tions for narrow band indices obtained by Worthey et al. 
(1994). In order to separate the age from metallicity we 
look at the properties of the galaxies in the four dimen- 
sional space of the parameters Fig, [MgFe], (1550 — V), 
and velocity dispersion E. The indices Hg and [MgFe] are 
particularly suited to separate age from metallicity effects 
because is a measure of the turn-off colour and lumi- 
nosity, and age in turn, whereas [MgFe] is more sensitive 
to the RGB colour and hence metallicity (Buzzoni et al. 
1992, 1993, Gonzales 1993). 

The paper is organized in the following way. Section 2 
briefly summarizes current information on basic data for 
elliptical galaxies, i.e. the observational hints on metal- 
licities and age ranges for the stellar content of ellipti- 
cal galaxies and bulges, the line strength indices Fig and 
[MgFe] the colour-magnitude relation (CMR) of elliptical 
galaxies and its current interpretation. Section 3 describes 
the recent models for elliptical galaxies by Bressan et al. 
(1994) and Tantalo et al. (1995). Section 4 presents the 
calculation of the narrow band indices for SSPs. Section 
5 deals with the integrated narrow band indices for mod- 
els of elliptical galaxies, with particular attention to the 
evolution in the Hg - [MgFe] plane. Section 6 addresses 
the question whether the observed distribution in the H^ 
- [MgFe] plane could result from recent bursts of star for- 
mation superposed to much older populations. Section 7 
analyses data and theoretical predictions in the four di- 
mensional space of coordinates H^, [MgFe], (1550 — V), 
and E, and presents a general discussion of the results 
emerging from this study. Section 8 deals with the gra- 



dients in Hj5 and [MgFe] across individual galaxies, and 
proposes a new method to separate age from metallicity. 
In particular it shows that different regions of a galaxy 
may have different age and metallicity and provides clues 
to understand the mechanism of galaxy formation. Finally, 
Section 9 presents some concluding remarks. 

2. Basic Observational data 

In this section we summarize current information on basic 
data for elliptical galaxies that are either constraints or 
targets of our models. 

2.1. Hints on metallicities and ages 

Direct determinations of metallicities in elliptical galax- 
ies are still rather limited. Most of the information comes 
from metal line-strength indices and their radial variations 
(Carollo et al. 1993, Carollo & Danziger 1994, Davies et 
al. 1993) or colour gradients (Schombert et al. 1993). The 
evidence arises that the metallicities for elliptical galax- 
ies are solar or larger. Furthermore, for bright elliptical 
galaxies there are also indications that the a-elements 
(O, Mg, Si, etc..) are enhanced with respect to Fe. The 
average [Mg/Fe], in particular, exceeds that of the most 
metal-rich stars in the solar vicinity by about 0.2-0.3 dex, 
and the ratio [Mg/Fe] is expected to increase with the 
galactic mass up to the this value. Due to their proxim- 
ity, the metallicity of stars in the Galactic Bulge is cur- 
rently measured by spectroscopic methods (cf. Rich 1990, 
McWilliam & Rich 1994). The mean metallicity is slightly 
lower than solar and long tails towards the metal-poor 
([Fe/H] = —1) and metal-rich end ([Fe/Hj = 1) are also 
present. High metallicities seem also to be indicated by 
the color-magnitude diagrams (CMD) of stellar popula- 
tions in the bulges of nearby Galaxies (M31 and M32). 
In the case of M32 (Freedman 1989, 1992), these stars 
are most likely in the RGB and AGB phases and metal- 
licities as high as [M/H]=0.1 are possible (cf. Bica et al. 
1990, 1991). Concerning the age, from the magnitude of 
the brightest AGB stars in M32, Freedman (1989, 1992) 
concludes that a fraction of these could be as young as 5 
Gyr. 

2.2. Line strength indices Fig and [MgFe] 

Table |l| contains the key data relevant to our purposes, 
namely the fully corrected indices Hp and [MgFe] and the 
velocity dispersion £ (km sec -1 ) of Gonzales (1993) for 
two different galaxy coverage, i.e. the central region within 
Re/ 8, where Re is the effective radius, and the wider region 
within Re/2. Column (1) is the galaxy identification NGC 
number, column (2) the logarithm of the Hg index, column 
(3) is the logarithm of the [MgFe] index, column (4) is the 
logarithm of the velocity dispersion £ in km sec -1 for the 
Re/8-data. Columns (5), (6) and (7) show the same but 
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Fig. 1. The fundamental plane Hp versus [MgFe] for the Re/8 
data of Gonzales (1993) together with the error bars for each 
galaxy 

for the Rc/2-data. The content of the remaining columns 
will be described below. 

The fundamental plane Hp versus [MgFe] for the Re/8 
data is shown in Fig. [l] together with the error bars for 
each galaxy. It is worth noticing that passing from the 
Re/8 to Re/2 data, there is a mean decrease in log [MgFe] 
amounting 0.032, and a mean increase in logflg amount- 
ing to 0.003. 

2.3. UV excess and (1550 - V) colour 

Much of the available information is from the study of 
Burstein et al. (1988). The main points are the following: 

(1) All studied elliptical galaxies have detectable flux 
short-ward of ~ 2000 A, with large galaxy to galaxy dif- 
ferences in the level of the UV flux. The intensity of the 
UV emission is measured by the colour (1550- V). 

(2) The colour (1550- V) correlates with the index Mg2, 
the velocity dispersion £ and the luminosity (mass) of the 
galaxy. The few galaxies (e.g. NGC 205) in which active 
star formation is seen do not obey these relations. 

(3) An important constraint is posed by the HUT observa- 
tions by Ferguson et al. (1991) and Ferguson & Davidsen 
(1993) of the UV excess in the bulge of M31. In this case 
the UV emission shows a drop-off short-ward of about 
1000 A whose interpretation requires that the tempera- 
ture of the emitting source must be about 25,000 K. Only 
a small percentage of the 912 < A < 1200 A flux can be 
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coming from stars hotter than 30,000 K and cooler than 
20,000 K. 

The (1550 - V) colours of Burstein et al. (1988) are 
listed in column (8) of Table [lj 

2.4- The CMR of early-type galaxies 

Colours and magnitudes of elliptical galaxies define a 
mean CMR, according to which the colours get redder at 
decreasing absolute magnitude and hence increasing mass 
of the galaxies. 

The slope of the CMR is commonly understood as 
indicative of a mass-metallicity sequence (Dressier 1984, 
Vader 1986) resulting from the effect of galactic winds 
(Larson 1974; Saito 1979a,b; Bressan et al. 1994). In 
brief, independently of their total mass, galaxies evolving 
as either closed or open boxes (cf. Tinsley 1980) should 
reach metallicities governed by the degree of gas gas con- 
sumption. Elliptical galaxies, having exhausted their gas, 
should possess almost identical (high) metallicities. There- 
fore, in order to have a mass-metallicity sequence, a mech- 
anism halting star formation and thus leaving different 
metallicities must be invoked: galactic winds are the right 
agent (cf. Arimoto & Yoshii 1986, 1987; Matteucci & Tor- 
nambe 1987; Yoshii & Arimoto 1987; Bressan et al. 1994; 
Tantalo et al. 1995). 

The tightness of the CMR seems to depend on the 
galaxy environment and could reflect the process of galaxy 
formation. 

Looking at the CMR of Bower et al. (1992a,b) for 
galaxies in the largest known nearby clusters Virgo and 
Coma, the colour dispersion is very small (uncertainty in 
absolute magnitudes is less of a problem here), with typi- 
cal rms scatter of 0.05 mag of which 0.03 mag can be ac- 
counted for by observational errors. Bower et al. (1992a, b) 
using the rate of colour change <9(U — Y)/dt for SSPs of 
Bruzual (1983) argued that for the Hubble time tu = 15 
Gyr and coordination in the galaxy formation process 
(their parameter (3 — 1) the galaxy ages are confined in 
the range 13.5 ~ tp ~ 15.0 Gyr. It is an easy matter to 
check that the Bower et al. (1992a, b) conclusion neither 
depends on the particular rate of colour change they have 
adopted nor the use of single SSPs instead of a manyfold 
of SSPs with different chemical composition weighted on 
the past rate of star formation (thus closer to the com- 
plexity of a real galaxy). This can be checked using the 
SSPs of Bressan et al. (1994) and Tantalo et al. (1995) or 
the models presented below. 

In contrast, examining the colour dispersion (relatively 
small also in this case) of the CMR for nearby galaxies in 
small groups and field, Schweizer & Seitzer (1992) found 
it compatible with recent mergers of spiral galaxies spread 
over several billion years. Also in this case the CMR im- 
plies a mass-metallicity sequence. 

Since the interpretation of the CMR bears very much 
on the kind of process of galaxy formation and evolution, 
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Table 1. Basic observational data used in this study: fully corrected indices H/3 and [MgFe], and velocity dispersion E (km 
sec -1 ) for the Re/8 and Re/2 data of Gonzales (1993); (1550 — V) colours of Burstein et al. (1988) for the central regions; 
magnitudes M v and M v , and colours (B-V) and (V-K) of Bender et al. (1992, 1993), Bower et al. (1992), Marsiaj (1992) and 
Schweizer & Seitzer (1992) as explained in the text. 
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its use as a constrain on the galactic models has different 
implications. If the CMR is indicating nearly coeval, old 
galaxies ranked along a mass-metallicity sequence, passive 
evolution after the initial star forming activity and onset 
of galactic winds is appropriate. This facilitates the use of 
the CMR as a constrain on galactic models (cf. Bressan 
et al. 1994, Tantalo et al. 1995). In contrast, if the CMR 
is compatible with mergers spreading over long periods of 
time, first the evolutionary history of an elliptical galaxy is 
more complicated, and the use of the CMR as a constraint 
is less secure. 

In order to cast light on this topic we try a closer in- 
spection of the data and derive the CMR for the Gonzales 
(1993) sample to be compared with that of Bower et al. 
(1992a, b) and Schweizer & Seitzer (1992). Cross-checking 
of the three lists reveals that, while the Gonzales (1993) 
and the Schweizer & Seitzer (1992) samples have many 
objects in common (17 galaxies), there is only one galaxy 



in common to all the three groups (NGC 4374), and four 
in common to Gonzales (1993) and Bower et al. (1992a, b), 
namely the galaxy above plus NGC 4472, NGC 4478, and 
NGC 4552. The galaxy NGC 4374 is particularly interest- 
ing because Schweizer & Seitzer (1992) have estimated for 
it a merger age of about 7 Gyr. 

In addition to this, we tried to assign magnitudes and 
colours to as many of the Gonzales (1993) galaxies as pos- 
sible making use of magnitudes and colours taken from 
Bender et al. (1992, 1993) and Marsiaj (1992). All the data 
are presented in Table [l]. Column (9) indicates whether 
the galaxy belong to the Schweizer & Seitzer (1992) list 
(SS) and/or the Virgo cluster (V or SSV); columns (10) 
and (11) are the Mb magnitude and the fine structure 
parameter S, respectively, of Schweizer & Seitzer (1992). 
The parameter S is a rough measure of dynamical youth 
or rejuvenation (0 no trace of fine structure, 10 strong fine 
structure). Column (12) and (13) are the magnitude Mb 
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and colour (B-V), respectively, taken from Bender et al. 
(1992, 1993). In a few cases these magnitudes are some- 
what different from those of Schweizer & Seitzer (1992). 
Columns (14) is the (B-V) colour listed by Marsiaj (1992). 
In general, it agrees with that of Bender et al. (1992, 1993). 
Columns (15) and (16) are the (V-K) colours of the Bower 
et al. (1992a, b) and Marsiaj (1992) lists, respectively. Also 
in this case, for the few galaxies in common the agreement 
is remarkable. This means that Marsiaj's (1992) (V-K) 
colours assigned to the remaining galaxies (when possi- 
ble) can be considered as reasonably good. Finally, col- 
umn (17) lists the magnitude My obtained from columns 
(12) and (13). 

The resulting (V — K) versus My CMR is shown in 
panel (a) of Fig. ||, whereas that of Bower et al. (1992a, b) 
is displayed in panel (b). In this latter, we adopt for Virgo 
the distance modulus (m — M)o=31.54 (Branch & Tam- 
mann 1992) and apply to Coma the shift <5(m — M)o = 
3.65 (cf. Bower et al. 1992a, b). In Fig. || we also plot the 
loci of constant ages (17, 10, 8 and 5 Gyr) for the models 
presented in Sect. 3. The comparison of the two CMRs 
yields the following results: 

(1) The slopes of the two CMRs agree each other and with 
the theoretical prediction. This means that both groups of 
galaxies reflect a similar mass-metallicity sequence. 

(2) Indeed the Bower et al. (1992a, b) CMR is less scattered 
than the one of Schweizer & Seitzer (1992), thus explaining 
the different interpretations. 

3) NGC 4374 and NGC 4697, for which rather recent 
merger and accompanying stellar activities have been pro- 
posed by Schweizer & Seitzer (1992), are compatible with 
the Bower et al. (1992a, b) conclusion at the same time. 

4) Finally, as galactic models are imposed to match the 
slope and the mean (V-K) colours of the CMR, the results 
do not depend on the particular CMR in usage. Therefore, 
owing to its better definition, we prefer to adopt the Bower 
et al. (1992) CMR to infer the slope of the underlying 
mass-metallicity sequence. 

3. Theoretical models of elliptical galaxies 

In this section we summarize the main properties of the 
models of elliptical galaxies used in the present analysis. 
They are taken from the studies of Bressan et al. (1994), 
Tantalo et al. (1995) in which an attempt is made to 
to explain in a coherent fashion the pattern of spectro- 
photometric properties of elliptical galaxies, i.e. the CMR, 
the origin of the UV excess and its dependence on the Mg2 
index and total luminosity, etc. The models in question 
differ for the kind of description adopted to follow the 
chemical history of the galaxy: the one-zone closed-box 
scheme in Bressan et al. (1994) and the one-zone infall 
scheme in Tantalo et al. (1995), this latter aimed at sim- 
ulating in a very simple fashion the formation of a galaxy 



Fig. 2. Panel a: the CMR for the sample by Schweizer & 
Seitzer for which we have obtained the magnitude My - The po- 
sition of NGC 4374 (filled star) and NGC 4697 (open square) 
is indicated. Panel b: the CMR for galaxies in Virgo (circles) 
and Coma (filled circles) reproduced from Bower et al.(1992). 
The galaxy NGC 4374 is shown by the open star. The typical 
error bars of the data are also displayed. In addition to this, 
in each panel, are drawn the line of constant age for the infall 
models characterized by the parameters: r = 0.10 Gyr, k = 1, 
£ = 0.50, and v variable with Ml (cf. Table 2). Each line corre- 
sponds to different values of the age: (solid: 17 Gyr), (dotted: 
10 Gyr), (short-dashed: 8 Gyr), and (long-dashed: 5 Gyr). See 
the text for more details 

by collapse of primordial gas (cf. Chiosi 1981 for more 
details). 

(1) Elliptical galaxies are supposed to be made of two 
components, i.e. luminous material of mass Ml embed- 
ded in a halo of dark matter of mass Md , whose presence 
affects only the gravitational potential of the galaxy and 
the binding energy of the gas. The key assumption of the 
closed-box models is that at the start of the star formation 
process (at time t=0) all the luminous mass is present in 
form of gas. In contrast, in the infall models the mass of 
the luminous component (supposedly in form of gas) is let 
increase with time according to 

^=M e XP (-t/r) (1) 

where r is the accretion time scale. The constant Mo is ob- 
tained from imposing that at the galaxy age Tq a certain 



(i 
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value of Ml(Tg) is reached. Therefore, the time depen- 
dence of Ml(£) is 



M L (t) 



M L (T G ) 



[1 - exp(-T G /r) 



■[l-exp(-t/r)] 



(2) 



As far as the dark component Md is concerned, we sim- 
ply assume it to be constant in time and equal to a cer- 
tain fraction of the asymptotic value of the luminous mass 
(Ml(Tq)/Md = /3), and to obey the spatial distribution 
with respect to Ml given by the model of Bertin et al. 
(1992) and Saglia et al. (1992), according to whom good 
fits of elliptical galaxies are possible for ratios Ml/Md 
and Rl/Rd equal to (3 = 0.2. With this assumption lit- 
tle effects are to be expected from the presence of dark 
matter. 

(2) The stellar birth rate. i.e. the number of stars of mass 
M born in the time interval dt and mass interval dM is 



dAf = tf(t, Z)(j){M)dtdM 



(3) 



where per ^(t) is the rate of star formation as a function 
of time and chemical enrichment, while <f>(M) is the initial 
mass function (IMF). The rate of star formation ^f(t) is 
taken to be proportional to the available gas mass M g (t) , 
i.e. 



= vMg(t) k M yr- 



(4) 



with k = 1. The constant of proportionality v represents 
the inverse star formation time scale in suitable units. The 
IMF is the Salpeter law 



4>{M) = M~ 



(5) 



with x — 2.35. The IMF is normalized by choosing the 
parameter Q 



J™" <p{M)xMxdM 



I M 



</>(M)xMxdM 



(6) 



fixing the fraction of total mass in the IMF above M* 
and deriving the lower limit of integration M^. The upper 
limit of integration is Mjj = 120M Q , the maximum mass 
in our data base of stellar models, while the mass limit 
M * is the minimum mass contributing to the chemical 
enrichment of the interstellar medium over a time scale 
comparable to the total lifetime of a galaxy. This mass is 
approximately equal to 1M Q . 

(3) The detailed treatment of the chemical enrichment, 
i.e. without instantaneous recycling (cf. Tinsley 1980), is 
included so that the temporal evolution of chemical abun- 
dances is properly taken into account (see Bressan et al. 
1994 for more details on the equations governing the chem- 
ical evolution of the models). 

(4) The models include the presence of galactic winds trig- 
gered by the energy deposit into the interstellar gas due 



Table 2. Key parameters of the infall models and correspon- 
dence between Ml and present-day total mass of stars Ms. 
Both are in units of 1O 12 M0. The time scale of mass accretion, 
t, and the age at which galactic winds occur, t gw , are expressed 
in Gyr. 



M L 


r 


V 


tgw 


Z 


< Z > 


M s 


C = 


0.50 












3.0 


0.10 


12.0 


0.26 


0.0710 


0.0360 


2.088 


1.0 


0.10 


7.2 


0.31 


0.0629 


0.0307 


0.648 


0.5 


0.10 


5.2 


0.35 


0.0543 


0.0265 


0.295 


0.1 


0.10 


3.0 


0.34 


0.0328 


0.0166 


0.041 


0.05 


0.10 


2.5 


0.29 


0.0235 


0.0122 


0.016 


0.01 


0.10 


1.0 


0.43 


0.0158 


0.0080 


0.002 



to the explosion of type I and II supernovae, and also stel- 
lar winds from massive stars. When the energy storage 
overwhelms the gravitational binding energy of gas this is 
supposed to be expelled from the galaxy thus halting fur- 
ther star formation. The prescription is exactly the same 
as in Bressan et al. (1994). 

(5) The models are identified by the mass of the luminous 
component Ml in units of 10 12 Mq (thereinafter referred 
to as M12). We remind the reader that this mass never 
coincides with the present day mass in stars. In the closed- 
box models Ml is the initial value of the luminous mass, 
part of which is converted into stars and part is expelled 
by galactic winds. Furthermore, the mass in living stars, 
Mg, most contributing to the light emitted by a galaxy is 
subjected to decrease with time because of the continuous 
death of stars leaving collapsed remnants and gas. In the 
case of infall models Ml is the asymptotic value of the 
luminous mass which is never reached because galactic 
wind is supposed to stop star formation and gas accretion 
at the same time (cf. Tantalo 1994, and Tantalo et al. 
1995). The same arguments for the mass in stars given for 
the closed box models apply also to the infall models. In 
Table || we give the correspondence between Ml and Mg 
at the present age (15 Gyr). 

(6) Finally, the isochrones at the base of the population 
synthesis calculations are from the Padua library (Bertclli 
et al. 1994). Ample ranges of chemical composition are 
considered, whose parameters Y and Z obey the enrich- 
ment law AY/AZ=2.5 (cf. Pagel 1989, Pagel et al. 1992). 
The main assumptions to be recalled here in view of the 
discussion below are (i) the use of the solar pattern of 
abundances for the so-called a-elements with respect to 
Fe, [a/Fe]=0; (ii) the use of rj = 0.45 in the mass loss rate 
along the RGB and AGB phases (cf. Bressan et al. 1994, 
Tantalo et al. 1995 for all details). 
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M L Age M bo i My (U - B) (B - V) (V - R) (V - K) (1550 - V) 



3.00 


17.0 


-24.514 


-23.662 


0.604 


0.999 


0.742 


3.330 


2.069 


3.00 


15.0 


-24.597 


-23.744 


0.597 


1.001 


0.743 


3.339 


2.217 


3.00 


12.0 


-24.766 


-23.916 


0.553 


0.981 


0.732 


3.347 


2.288 


3.00 


10.0 


-24.972 


-24.110 


0.525 


0.968 


0.727 


3.375 


2.155 


3.00 


8.0 


-25.146 


-24.315 


0.483 


0.937 


0.711 


3.339 


2.632 


3.00 


5.0 


-25.560 


-24.741 


0.444 


0.911 


0.691 


3.348 


6.093 


1.00 


17.0 


-23.282 


-22.470 


0.594 


1.000 


0.735 


3.273 


2.715 


1.00 


15.0 


-23.369 


-22.556 


0.586 


1.002 


0.736 


3.280 


2.938 


1.00 


12.0 


-23.542 


-22.733 


0.547 


0.986 


0.726 


3.289 


3.160 


1.00 


10.0 


-23.738 


-22.919 


0.522 


0.973 


0.720 


3.316 


3.070 


1.00 


8.0 


-23.918 


-23.122 


0.479 


0.946 


0.705 


3.287 


3.536 


1.00 


5.0 


-24.332 


-23.551 


0.420 


0.908 


0.682 


3.280 


6.068 


0.50 


17.0 


-22.448 


-21.668 


0.573 


0.994 


0.729 


3.222 


3.396 


0.50 


15.0 


-22.537 


-21.758 


0.566 


0.997 


0.729 


3.225 


3.792 


0.50 


12.0 


-22.714 


-21.940 


0.531 


0.983 


0.720 


3.234 


4.595 


0.50 


10.0 


-22.901 


-22.118 


0.507 


0.971 


0.713 


3.258 


4.712 


0.50 


8.0 


-23.084 


-22.320 


0.465 


0.946 


0.699 


3.233 


5.096 


0.50 


5.0 


-23.498 


-22.751 


0.396 


0.902 


0.674 


3.217 


6.019 


0.10 


17.0 


-20.362 


-19.652 


0.472 


0.947 


0.706 


3.080 


3.359 


0.10 


15.0 


-20.454 


-19.752 


0.470 


0.955 


0.706 


3.071 


3.887 


0.10 


12.0 


-20.639 


-19.945 


0.445 


0.946 


0.698 


3.073 


5.402 


0.10 


10.0 


-20.802 


-20.110 


0.420 


0.931 


0.689 


3.078 


5.937 


0.10 


8.0 


-20.994 


-20.314 


0.384 


0.909 


0.677 


3.064 


6.061 


0.10 


5.0 


-21.398 


-20.730 


0.321 


0.862 


0.652 


3.051 


5.804 


0.05 


17.0 


-19.337 


-18.675 


0.416 


0.915 


0.691 


2.977 


3.243 


0.05 


15.0 


-19.433 


-18.780 


0.418 


0.928 


0.692 


2.969 


3.818 


0.05 


12.0 


-19.620 


-18.974 


0.397 


0.922 


0.686 


2.970 


5.373 


0.05 


10.0 


-19.778 


-19.137 


0.373 


0.908 


0.677 


2.970 


5.915 


0.05 


8.0 


-19.979 


-19.344 


0.342 


0.887 


0.666 


2.963 


5.993 


0.05 


5.0 


-20.359 


-19.741 


0.279 


0.836 


0.639 


2.941 


5.628 


0.01 


17.0 


-17.282 


-16.704 


0.329 


0.860 


0.664 


2.789 


3.037 


0.01 


15.0 


-17.393 


-16.819 


0.335 


0.880 


0.667 


2.793 


3.751 


0.01 


12.0 


-17.580 


-17.014 


0.317 


0.879 


0.664 


2.794 


5.381 


0.01 


10.0 


-17.740 


-17.176 


0.298 


0.867 


0.656 


2.798 


5.805 


0.01 


8.0 


-17.957 


-17.393 


0.275 


0.850 


0.647 


2.803 


5.792 


0.01 


5.0 


-18.307 


-17.768 


0.214 


0.792 


0.616 


2.757 


5.225 



3.1. Main results for infall models 

For the reasons amply discussed by Bressan et al. (1994) 
and Tantalo et al. (1995) infall models match the region 
2000 - 3500 A of the spectrum of an elliptical galaxy much 
better than the closed-box ones. Therefore in the anal- 
ysis below we adopt the infall models of Tantalo et al. 
(1995). They successfully reproduce the CMR for Virgo 
and Coma galaxies by Bower et al. (1992a,b) and the UV 
properties of elliptical galaxies for the following parame- 
ters: enrichment law AY/AZ=2.5, infall time scale r = 0.1 



Gyr, k = 1, £ = 0.50, and v decreasing from 12 for the 3 
M12 galaxy down to 2.5 for the 0.05 M12 object. 

The main properties of these models are presented in 
Table ^ which contains the specific star formation effi- 
ciency the age at the onset of the galactic wind, the 
maximum metallicity Z, the mean metallicity < Z >, and 
the mass Ms in stars at the present epoch. 

The chemical structure of the models is represented 
by the normalized, cumulative distribution of the mass 
in stars per metallicity bin at at the present epoch (say 
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15 Gyr) displayed in Fig. [| It is worth noticing the very 
small percentage of stars of low metallicity as compared to 
the case of the closed-box models of Bressan et al. (1994). 



Fig. 3. Chemical structure of the galactic models with infall 
and different Ml as indicated. The curves are the cumulative 
fractionary mass of living stars per metallicity bin normalized 
to the present-day value of mass in stars (Ms). The parameters 
of the models are those of Table 2. See the text for more details 

A summary of the temporal evolution of absolute 
bolometric and visual magnitudes, Mboi and My respec- 
tively, the broad-band colours (U - B), (B - V), (V - R), 
(V — K), and colour excess (1550 — V) is given in Table | 
for selected values of the age. The CMR of these models is 
shown in Fig. || for various values of the age as indicated. 

Like in Bressan et al. (1994) UV excess is generated by 
a suitable admixture of three components (see also Greg- 
gio & Renzini 1990). 

(1) The classical post asymptotic giant branch (P-AGB) 
stars (see Bruzual 1992, Bruzual & Chariot 1993, Chariot 
& Bruzual 1991) P-AGB stars are always present in the 
stellar mix of a galaxy. The major problem with these 
stars is their high effective temperature and the relation 
between their mass and that of the progenitor. The initial- 
final mass relationship is not firmly established. The most 
popular empirical determination is by Weidemann (1987). 
Matching this relation depends on the efficiency of mass 
loss during the AGB phase and other details of model 
structure (cf. Chiosi et al. 1992). 

(2) The hot horizontal branch (H-HB) and AGB-manque 
stars of very high metallicity (say Z ~ 0.07). which are 
expected to be present albeit in small percentages in the 



stellar content of bulges and elliptical galaxies in general. 
Indeed, these stars have effective temperatures in the right 
interval and generate ISEDs whose intensity drops short- 
ward of about 1000 A by the amount indicated by the ob- 
servational data. (Ferguson et al. 1991, Ferguson & David- 
sen 1993). The formation of these stars can occur either 
with low values of the enrichment ratio {AY/ AZ ~ 1) and 
strong dependences of the mass-loss rates on the metallic- 
ity (Greggio & Renzini 1990) or even with canonical mass- 
loss rates and suitable enrichment laws {AY/ AZ ~ 2.5) 
as in Horch et al. (1992) and Fagotto et al. (1994a, b,c). 

(3) Finally, the very blue HB stars of extremely low metal- 
licity (Lee 1994). These stars have effective temperatures 
hotter than about 15,000 K but much cooler than those 
of the P-AGB stars. Therefore, depending on their actual 
effective temperature, they can generate ISEDs in agree- 
ment with the observational data. However, most likely 
they are not the dominant source of the UV flux because 
the analysis by Bressan et al. (1994) clarifies that in the 
wavelength interval 2000 < A < 3000 A the ISEDs of 
the bulge of M31 and of elliptical galaxies like NGC 4649 
are fully consistent with the notion that virtually no stars 
with metallicity lower than Z = 0.008 ought to exist in 
the mix of these stellar populations. 

Our model galaxies with AY/AZ=2.5 emit UV radia- 
tion by old H-HB stars only for ages older than about 5.6 
Gyr and in presence of even a tiny fraction of stars in suit- 
able metallicity bins (say Z > 0.07). The precise value of 
the age depends also on the final-initial mass relationship 
for the emitting stars (cf. Bressan et al. 1994, Tantalo et al. 
1995 for details). Fig. || show the (1550 - V) colour of our 
galactic models as a function of the age. It is worth point- 
ing out the different behaviour of the (1550 — V) colour 
passing from the 3 M12 galaxy, in which H-HB, AGB- 
manque, and P-AGB stars concur to generate the UV 
flux, to the 0.01 M12 galaxy, in which only P-AGB stars 
produce the UV radiation. This trend is caused by the 
decreasing mean and maximum metallicity at decreasing 
galaxy mass (cf. Bressan et al. 1994 for details). 

Concluding this section, it is worth pointing out that in 
our models the efficiency of star formation per unit mass 
of gas increases at increasing galactic mass {v goes from 
2.5 to 12.0 as the mass increases from 0.05 to 3 M12). This 
is the reason for the nearly constant value of age t gw at 
which galactic winds occur halting star formation. Similar 
trend for the parameter v has been invoked by Matteucci 
(1994) to explain the abundance ratios [Mg/Fe] observed 
in elliptical galaxies (see below). 

3.2. General remarks on the models 

We like to call the attention on three somewhat critical 
aspects of the models that could be used to invalidate the 
present analysis. The points in question are (i) the adop- 
tion of the one-zone description which does not allow us to 
take into account spatial gradients in colours and metal- 
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Age (Gyr) 

Fig. 4. The time variation of the (1550 — V) colour of the 
galactic models with infall and different Ml as indicated. The 
parameters of the models are those of Table 2. See the text for 
more details 

licity, (ii) the solar partition of abundances used in the 
stellar models, and finally (iii) our treatment of galactic 
winds. 

(i) It is clear that the (V — K) colours and magni- 
tudes of the galaxies in the Bower et al. (1992a, b) sample 
refers to the whole galaxy, whereas the UV excess and its 
companion colour (1550 — V) of Burstein et al. (1988) in 
most cases refer to the central region of a galaxy. This 
is also true for the various indices measured by Gonzales 
(1993) to be discussed below. The existence of gradients in 
colours and metallicities across elliptical galaxies is a well 
established observational fact (Carollo et al. 1993; Car- 
olio & Danziger 1994, Davies et al. 1993, Schombert et al. 
1993). This means that insisting on the one-zone descrip- 
tion as in the present models, could be source of difficulty 
when comparing model results with observational data. It 
seems reasonable to argue that higher metallicities can be 
present in the central regions thus facilitating the forma- 
tion of the right type of stellar sources of the UV radiation 
(H-HB and AGB manque stars), whereas lower metallici- 
ties across the remaining parts of the galaxies would make 
the other integrated colours such as (V-K) bluer than ex- 
pected from straight use of the one-zone model for the 
whole galaxy. 

(ii) Recent observations indicate that the pattern of 
abundances in elliptical galaxies is skewed towards an 
overabundance of a-elements with respect to Fe (cf. Car- 
ollo et al. 1993; Carollo & Danziger 1994). As already 
recalled, the library of stellar models in use is based 



on the standard (solar) pattern of abundances. Work is 
in progress to generate libraries of stellar models and 
isochrones with [a/Fe] > 0. Preliminary calculations for 
the solar metallicity (Z=0.02) show the effect to be small 
in the theoretical CMD. In addition, calculations of the 
Mg2 index from synthetic spectra for different metallici- 
ties and partitions of a-elements (Barbuy 1994) show that 
passing from [a/Fe = 0.0 to [cn/Fe = 0.3 the Mg2 index of 
an old SSP (15 Gyr) increases from 0.24 to 0.30. 

(iii) Concerning galactic winds, we would like to com- 
ment on the recent claim by Gibson (1994) that Bressan 
et al. (1994) underestimate the effect of supernova explo- 
sions and on the contrary overestimate the effect of stellar 
winds from massive stars. In brief, Bressan et al. (1994) as- 
suming the CMR to be a mass-metallicity sequence looked 
at the metallicity that would generate the right colours. 
They found that considering supernovae as the only source 
of energy, by the time the galactic winds occur, the gas 
fraction has become too low and the metallicity too high 
in turn, that the CMR is destroyed (it runs flat). To cope 
with this difficulty they included another source of energy, 
i.e. the stellar winds from massive stars. In analogy with 
the behaviour of a SN remnant's thermal energy, which 
decreases with time as e$N °c (t/t c )~°- 62 for times greater 
than t c , Bressan et al. (1994) assumed that the energy 
of stellar winds obeys the same law but with a different 
t c , i.e. e oc (t/t cw )~ - 62 for times greater than t cw . The 
time t cw was considered as a parameter. They found that 
t cw = 1.5 x 10 7 yr is a good choice. This time roughly 
corresponds to the evolutionary lifetime of a 10 M Q star. 
In other words it is the time scale over which a group of 
newly formed O-type stars would evolve away from the 
SSP. They found that t cw shorter than this would not 
allow sufficient powering of the interstellar medium, so 
that galactic winds would occur much later than required 
by the CMR. Bressan et al. (1994) provided also an ex- 
planation for the different results found by other authors 
(Arimoto & Yoshii 1987, Matteucci & Tornambc' 1987, 
Angeletti & Giannone 1990, and Padovani & Matteucci 
1993). Admittedly, this additional source of energy was 
invoked to keep the standard interpretation of the CMR. 

4. Narrow-band indices for SSPs and galaxies 

In this paper we make use of the recent library of empiri- 
cal calibrations of twenty one indices of absorption feature 
strengths calculated by Worthey (1992) and Worthey et 
al. (1994). The calibrations are presented in form of an- 
alytical fits that give the index strength as a function of 
stellar atmospheric parameters, i.e. e = 5040/T e ff, \ogg 
(gravity), and [Fe/H] both for warm (up to 13,260 K) and 
cool stars (down to 3570 K). The definition of the indices 
are given in Worthey (1992) and Worthey et al. (1994) 
to whom we refer. Suffice the recall here that the index 
[MgFe] is the geometric mean of the two indices Mgb and 
< Fe >, [MgFe] = (Mgb < Fe >) 5 where in turn the in- 
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dcx < Fc > is the mean of the indices Fes27o and Fes335, 
< Fe > = (Fe 52 7o + Fe 533 5)/2. 

According to their definition, all indices are con- 
structed by means of a central band-pass and two pseudo- 
continuum band-passes on either side of the central band 
(cf. Worthey et al. 1994 for details). The continuum flux 
is interpolated between the mid points of the pseudo- 
continuum band-passes. 

The integrated indices for SSPs and model galaxies are 
calculated using the following method. 

For every combination of T e g, g and [Fe/H], first we 
derive the flux in the continuum band-pass from the li- 
brary of stellar spectra used by Bressan et al. (1994) im- 
plemented by the revision of Tantalo et al. (1995), and 
then we calculate the flux in the central pass-bands with 
the analytical fits of Worthey et al. (1994) for the same 
values T e ff, g and [Fe/H]. Let us call Fj\ and Fcx the flux 
in central pass-band and in the continuum, respectively. 

The integrated Fix and Fcx (thereinafter indicated as 
F K ) for the stellar content of a galaxy of age T are given 

by 

F K {T)= / V(t,Z)cj>{M) F K (M,r',Z)dtdM (7) 
Jo J M L 

where F k (M,t', Z) is one of the two fluxes for a star of 
mass M, metallicity Z(t), and age t' = T — t. Separat- 
ing the contribution from single SSPs, the above integral 
becomes 

F K {T) = [ q/(t,Z) F K . ssp (r\Z) dt (8) 
Jo 

where 

r-Mu 

F KiSSP (t',Z)= <j>{M) F k (M,t',Z) dM (9) 

J Ml 

are defined as the integrated fluxes, either Fix or Fcx 
for a single SSP, i.e. of a coeval, chemically homogenous 
assembly of stars with age r' and metallicity Z. 

Knowing the integrated fluxes either for a SSP or a 
galaxy, the definition of each index is applied to get back 
the integrated index. 

The rate of star formation ^(t, Z) and initial mass 
function 4>(M) in equation (|^) are the same as in equations 
(§, (|), and (§). 

Tabulations of indices for SSPs with different chemi- 
cal composition and age are given in Table 4. In order to 
present results independent of the particular choice made 
for the IMF normalization, the SSPs of Table 4 are calcu- 
lated assuming the Salpeter IMF with x = 2.35, and fixed 
lower and upper limits of integration, Ml = 0.15Mq and 
Mu = 120M Q , respectively. Care must be paid applying 
these indices for SSPs to real assemblies of stars (either 
clusters or galaxies). 

The correlation between tig and [MgFe] for SSPs with 
different metallicities is shown in Fig. ||. The age goes from 




Fig. 5. Single Stellar Populations of different metallicity in the 
H/3 - [MgFe] plane. Along each the age varies from 1 to above 
15 Gyr as indicated. The full dots and open triangles are the 
Re/8-data and Re/2-data of Gonzales (1993). The error bars 
of the data are those of Fig. ^ 

1 Gyr to beyond 15 Gyr in steps of A logt = 0.05 as indi- 
cated. The inspection of the data in Table 4 and compari- 
son with similar results obtained by Gonzales (1993) show 
that excellent agreement exists between the two studies. 
The dip and the loop-like structure in the SSPs with the 
lowest and highest metallicity are caused by the onset of 
the HB and appearance of the H-HB stars, respectively. 
Only for the sake of provisional comparison we plot in 
Fig. U also the Re/2-data (open triangles) and Re/8-data 
(filled circles) of the Gonzales (1993) sample. It is soon 
evident that the locus drawn by the observations has the 
same slope of the SSPs with metallicity comprised between 
Z=0.02 and Z=0.05. 

5. Elliptical galaxies in the - [MgFe] plane 

With the aid of the technique described in the previous 
sections, we calculate the temporal evolution of the nar- 
row band indices for model galaxies with different Ml- 
Table || lists the indices rig [MgFe] and Mg2 of our galac- 
tic models at selected values of the age, namely 17, 15, 12, 
10, 8 and 5 Gyr. 

The evolutionary path of model galaxies in the EL3 - 
[MgFe] plane is shown in Fig. |] together with the Re/8- 
data of the Gonzales (1993) sample. The galactic models 
are labelled by their Ml . Along each sequence the age in- 
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M L 


A PY 1 






loeMe-2 


M L 


Ase 


lop" l~T/a 


loe- [MeFel 


loe-Me-2 


3.00 


17.0 


0.157 


0.552 


0.277 


3.00 


10.0 


0.204 


0.533 


0.261 


1.00 


17.0 


0.167 


0.536 


0.265 


1.00 


10.0 


0.209 


0.517 


0.250 


0.50 


17.0 


0.177 


0.520 


0.254 


0.50 


10.0 


0.215 


0.500 


0.239 


0.10 


17.0 


0.211 


0.468 


0.221 


0.10 


10.0 


0.241 


0.446 


0.205 


0.05 


17.0 


0.234 


0.433 


0.201 


0.05 


10.0 


0.256 


0.412 


0.186 


0.01 


17.0 


0.278 


0.368 


0.170 


0.01 


10.0 


0.284 


0.349 


0.157 


3.00 


15.0 


0.155 


0.551 


0.277 


3.00 


8.0 


0.245 


0.518 


0.249 


1.00 


15.0 


0.163 


0.535 


0.265 


1.00 


8.0 


0.245 


0.502 


0.238 


0.50 


15.0 


0.172 


0.518 


0.253 


0.50 


8.0 


0.248 


0.485 


0.228 


0.10 


15.0 


0.201 


0.466 


0.219 


0.10 


8.0 


0.273 


0.430 


0.195 


0.05 


15.0 


0.220 


0.432 


0.200 


0.05 


8.0 


0.289 


0.396 


0.177 


0.01 


15.0 


0.257 


0.369 


0.169 


0.01 


8.0 


0.315 


0.334 


0.150 


3.00 


12.0 


0.179 


0.539 


0.267 


3.00 


5.0 


0.298 


0.496 


0.232 


1.00 


12.0 


0.185 


0.524 


0.255 


1.00 


5.0 


0.308 


0.476 


0.220 


0.50 


12.0 


0.191 


0.508 


0.244 


0.50 


5.0 


0.317 


0.457 


0.209 


0.10 


12.0 


0.216 


0.457 


0.212 


0.10 


5.0 


0.345 


0.400 


0.179 


0.05 


12.0 


0.232 


0.423 


0.193 


0.05 


5.0 


0.362 


0.362 


0.161 


0.01 


12.0 


0.261 


0.360 


0.163 


0.01 


5.0 


0.394 


0.291 


0.133 



creases from the top to the bottom, going from 1 to beyond 
15 Gyr as indicated. Although the Re/2-data would indeed 
be closer to the theoretical results for the one-zone models 
we are using, we prefer to adopt the Re/8-data referring 
to the central part of the galaxies. However, we would like 
to clarify that similar results would be obtained using the 
Re/2-data. Indeed passing from Re/2 to Re/8 the net ef- 
fect is that 13,3 and [MgFe] indices get slightly " redder" by 
the quantities A log fig = 0.003 and A log [MgFe] = 0.032 
as indicated by the vectors shown in Fig ||. 

In any case, compared with the Re/8-data, our values 
for [MgFe] are somewhat bluer than observed. This could 
reflect the solar partition of abundances adopted for the 
stellar models in use. It could also mean that the very 
central parts of the galaxies are more metal-rich than pre- 
dicted by the one-zone model. Indeed there is much better 
agreement between the Re/2-data and theoretical mod- 
els. Therefore, we consider the above discrepancy to be 
marginal and not invalidating the analysis below. It will 
be taken into account whenever necessary. On the basis of 
the above remarks, we consider it reasonable to impose co- 
incidence between our reddest models and the red edge of 
the observational distribution, and apply to the models an 
offset of A log [MgFe] = 0.05. The required shift is compa- 
rable with the increase in the Mg2 index found by Barbuy 
(1994) for [Mg/Fe] = 0.3 and reported in section 3.2. This 
shift will always be applied when comparing models with 
observations. The same is true for the Re/2-data. In such 
a case the shift along the [MgFe] axis amounts only to 
about A log [MgFe] = 0.02. 



As expected the evolutionary paths of the galaxies run 
parallel to that of SSPs. This is the consequence of in- 
terpreting the CMR as a mass-metallicity sequence, along 
which both the mean and the maximum metallicity of a 
galaxy are ultimately determined by the onset of galactic 
winds. In the canonical view of the CMR (slope), galac- 
tic winds and consequent interruption of the star forming 
activity occur earlier at decreasing galactic mass (lower 
mean and maximum metallicities). Since the subsequent 
evolution is passive at frozen chemical composition, the 
behaviour of a SSP is mimicked. Therefore, in the Ebj 
- [MgFe] plane galactic models shift to the left at de- 
creasing mass and mean metallicity in turn, cf. the en- 
tries of Table || and the values marked along the line of 
Fig. ^ labelled " CMR-strip" (see below) which shows the 
mean slope of the old isochrones (10 to 17 Gyr). The 
mean metallicity is annotated above, whereas the maxi- 
mum metallicity is written below the line in correspon- 
dence to the galactic models with different Ml- 

For the purposes of comparison, in Fig. ^| we also 
show the most massive closed-box model of Brcssan et 
al. (1994), namely the case with 3 M12 (triangles). It is 
interesting to note the large effect passing from closed- 
box to infall models. At given Ml, the closed-box models 
are much bluer both in tig and [MgFe] than the infall 
models. No galaxy is found in the region occupied by the 
closed-box model with 3 M12 and lower. This is another 
argument in favour of the infall models. 

In order to quantify the effect of variations in metallic- 
ity and age on the position of a galaxy in the EL3 - [MgFe] 
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plane, we calculate the vectors A log Hp and A log [MgFe] 
at constant age and constant metallicity. To this aim, we 
consider models of given Ml to derive (AlogH / g)< and 
(A[MgFe])f and models of fixed age and different Ml to 
get (AlogH^)^ and (A log [MgFe] )z, with obvious mean- 
ing of the symbols. We find the following relations over 
wide ranges of ages and metallicities 



(AlogH,) t = (^) t .-1.7 



(A log H-p)z 
and 



, d log H 



0- 



dt 



-0.01 Gyr' 1 



(Alog[MgFe]) t = (^i»l jt 



(A log [MgFe]) z 



dZ 



.dlogjMgFe] 

( — at — )z 



3.5 



(10) 



(11) 



(12) 



0.005 Gyr' 1 (13) 



Looking at the data of Fig. |(| although some scatter 
is present along the [MgFe] axis, most of the galaxies fall 
into the region between the sequences with 0.5 and 3 M 12 . 
This implies that the metallicity distribution within most 
of the galaxies in the sample is quite similar with both 
the mean and maximum metallicity nearly identical (see 
the metallicity distributions of the 0.5 and Ml galaxies 
shown in Fig. |3|). However, there are a few galaxies (M32 
and NGC 7454, in particular) that apparently fall along 
sequences of smaller Ml and hence, in our scheme, lower 
metallicities. In contrast, there is a large scatter along the 
Hp axis that seems to hint large variations in the age from 
galaxy to galaxy. This indeed is the conclusion reached by 
Gonzales (1993). For instance, the prototype galaxies M32 
and NGC 4649 would turn out to be as young as 3 Gyr 
and as old as 17 Gyr, respectively. 

There is a puzzling feature coming out from the close 
inspection of Fig. 0. In brief, as the galactic models have 
been constrained to match the slope and the (V — K) 
colours of the CMR under the assumption that this latter 
a metallicity-mass sequence of old (say 15 Gyr), nearly 
coeval objects, we expect galaxies to be located in the 
Hp - [MgFe] plane along a the line indicated as CMR- 
strip, which is the analog of the CMR of Fig 0. It is soon 
evident that the slope of the CMR-strip does not agree 
with the slope traced by the observational data. One may 
argue that this results from the lack of homogeneity be- 
tween the samples of Gonzales (1993) and Bower et al. 
(1992a, b). Coevality holds for the cluster galaxies of Bower 
et al. (1992a, b), whereas the spread in age applies to the 
field galaxies of Schweizer & Seitzer (1992) and Gonzales 
(1993). In both cases, as the CMRs of cluster and field 
galaxies have similar slope, one would expect galaxies to 
scatter within the area comprised between the youngest 
and oldest isochrone pertinent to the sample. 



Q5 1 0.5 1 3 



R e /8 aperture 



a LogH, 
a l,og[MgFe] 



same aae 
same mass 



CMR strip 
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Fig. 6. Evolution of the infall models in the Hp - [MgFe] 
plane (solid lines). Each model is labelled by Ml- The dot- 
ted lines are isochrones of different age (in Gyr) as indicated. 
The line called CMR-strip with the mean and maximum metal- 
licities of the galactic models annotated along it is the analog 
of the CMR. The solid triangles show the evolutionary path 
of the 3 M12 galaxy of Bressan et al. (1994) calculated with 
the closed-box formalism. The full dots are the Re/8-data of 
Gonzales (1993) . Finally, the vectors A log H/j and A log [MgFe] 
display the mean shift in the data going from the R/8-data to 
the Re/2-data. The error bars of the data are those of Fig. jj] 
See the text for more details 



In the scenario of coeval, old galaxies, the four objects 
of Bower et al. (1992a, b) in common with Gonzales (1993) 
- see the entries of Table || - should fall in the red corner 
of the H<3 - [MgFe] plane with little dispersion. This is 
not true, because NGC 4478 has a significantly bluer Hp 
suggesting a younger age, much younger than the range 
permitted by the tightness of the CMR. Conversely, in 
the scenario of galaxies with different ages (as perhaps 
induced by mergers), the associated chemical enrichment 
must be such that galaxies acquire very similar metallicity 
distributions. 

What we learn from this preliminary analysis of the 
observed distribution of galaxies in the Hp - [MgFe] plane 



In spite of their total luminosity and mass, most galax- 
ies seem to possess nearly identical chemical struc- 
tures. 

Ages seem to vary from galaxy to galaxy. 
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— Galaxies do not simply distribute along the locus ex- 
pected for objects matching the mass-metallicity se- 
quence of the CMR. 

The following questions can be addressed: Is the age 
spread real ? Does it imply a truly different starting epoch 
of the star formation activity, or other effects could change 
the evolutionary path of a galaxy in the Hg - [MgFe] plane 
and mimic a spread in age ? Why the predictions from the 
CMR do not find close correspondence in the H^ - [MgFe] 
plane ? 



• » 




t 



Fig. 7. Effects of a recent burst of star formation on the evolu- 
tionary path in the H/3 - [MgFe] plane. The galaxy has mass of 
3 Mi 2. The burst is supposed to start at 12 Gyr with such an 
intensity that 1.0% of the luminous mass of the galaxy in form 
of gas is turned into stars. The full dots are the Re/8-data of 
Gonzales (1993). The error bars of the data are those of Fig. nl 



6. Recent episodes of star formation ? 

As already recalled, a growing body of literature on fine 
structures and kinematic anomalies suggest the early type 
galaxies were either formed or restructured by mergers 
in relatively recent times (cf. Schweizer & Seitzer 1992, 
Chariot & Silk 1994, and references). In addition, studies 
of population synthesis (Faber et al. 1992 and references) 
and photometric data (color-magnitude diagrams) of the 
stellar content in nearby galaxies (cf. the case of M32 by 
Freedman 1989, 1992) seem to suggest that E+SO galax- 
ies contain various admixtures of intermediate age stars. 



Fig. 8. The variation of the H,3 index as a function of time 
for the galaxy with a recent burst of star formation shown in 
Fig. [?] The thin lines show two SSPs with different metallicities 
as indicated 

The merger hypothesis implies star burst of various in- 
tensity (cf. the models of merging spirals by Alfensleben 
& Gerhard 1994), whose net result is the formation of 
an E/SO galaxy. It is worth clarifying here that the term 
"merger" may refer to quite different views: either the 
classical merger of equal-mass spiral galaxies (Schweizer 
& Seitzer 1992) or most of the stars forming at an early 
epoch and very modest star formation at later epochs (low 
red-shifts) as pointed out by Chariot & Silk (1994), who 
apparently concluded that the data could not help to dis- 
tinguish between these two scenarios. 

How a recent burst of star formation, superposed to 
much older populations, would affect the evolutionary 
path of a galaxy in the Hg - [MgFe] plane is easy to foresee 
by means of the following experiments. Having in mind the 
case of M32 for which the existence AGB stars as young 
as about 3 -j- 5 Gyr (Freedman 1989, 1992) has been sug- 
gested, we suppose that at the arbitrary age of 12 Gyr a 
small burst of short duration occurs. The intensity of the 
burst is such that only 1.0% of the current luminous mass 
of the galaxy is turned into stars, whereas its duration 
is only 1 x 10 8 yr. As expected and in analogy to what 
is known from similar experiments in the (U — B) and 
(B — V) colours (cf. the models of Alfensleben & Gerhard 
1994), the occurrence of a burst of star formation induces 
a loop in the H3 - [MgFe] plane. At the start of the burst 
both H^ and [MgFe] get bluer, and when the burst is over 
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they recover the previous values. This is shown in Fig. [7] 
and Fig. ||. For this burst, the recovery time scale is about 
1 Gyr. Therefore, in order to be detected a burst of this 
type should have occurred as recently as 1 Gyr ago. Of 
course, the intensity and duration of the burst may change 
from galaxy to galaxy, so that both the amplitude of the 
loop in the - [MgFe] plane and the time scale at which 
the pre-burst values are recovered may be different from 
the ones we have been using. Stronger and longer bursts 
induces larger amplitudes of the loop and somewhat longer 
recovery time scales. 

On the basis of these experiments, one could argue 
that the large spread along the Hp direction is the result 
of a recent episode of star formation (either triggered by a 
merger or caused by an internal agent) superposed to an 
older population whose properties are temporarily masked 
by the star forming activity so that deciphering the true 
age of the system from the location in the EFj - [MgFe] 
plane is not possible. 

However, before accepting this conclusion, we feel it 
wise to examine in detail the correlation of and [MgFe] 
with the colour (1550 — V) and the velocity dispersion S. 




Fig. 9. The relation between (1550 — V) and The indices 
are from Gonzales (1993), whereas the colours (1550 — V) 
are from Burstein et al. (1988). The crosses show the uncer- 
tainty affecting the observational data. The solid and dotted 
lines correspond to loci of constant age (in Gyr) and mass (in 
units of M L as indicated. The position of M32 and NGC 584 
is brought into evidence 




Fig. 10. The relation between (1550 - V) and [MgFe]. The 
indices [MgFe] are from Gonzales (1993), whereas the colours 
(1550 — V) are from Burstein et al. (1988). The various sym- 
bols have the same meaning as in Fig. ^| 

7. The four-dimensional space - [MgFe] - 

(1550 - V) and £ 

The narrow range of the mean and maximum metallici- 
ties inferred from the - [MgFe] plane (with the few ex- 
ceptions to be discussed separately), their relatively high 
values suggested by the comparison with theoretical mod- 
els (see also the summary on observational hints), and 
the provisional hypothesis that all galaxies are old nearly 
coeval objects would imply that in contrast to the ob- 
servations (see Table 0) all galaxies of our sample should 
emit UV radiation at comparable levels (see the variation 
of the (1550 — V) colour as a function of the age for the 
0.5 and 3 Ml galaxies). Unless galaxies are younger than 
a few Gyr, this remark would apply also to the case in 
which a significant age difference from galaxy to galaxy is 
allowed. 

The situation is shown in the Fig. || and [lO] display- 
ing the separated relations Ebj - (1550 — V) and [MgFe] 
- (1550 — V) for the galaxies in Table |l|. As well known 
the (1550 — V) colour gets bluer (stronger UV flux) with 
increasing [MgFe] and decreasing Hg. The data are com- 
pared with the locus of theoretical models with different 
Ml, and mean metallicity in turn, and different age. The 
solid and dashed lines in Figs. || and [To] show the loci of 
constant mass and constant age, respectively. The compar- 
ison is made using the Re/8-data because the (1550 — V) 
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Table 6. Mass in stars, Ms, radius of the luminous component Rl, gravitational potential Q, and velocity dispersion E of the 
model galaxies. The units of each quantity are indicated. 



M L 


M s 


Rl 


Q 


E 


log'S 


log E ad 


W^Mq 


10^M© 


Kpc 


ergs 


Kmsec~ L 






3 


2.08 


39.1 


6.18(60) 


308 


2.488 


2.506 






47.8 


1.05(61) 


334 


2.523 




1 


0.65 


20.6 


1.14(60) 


237 


2.374 


2.395 






26.1 


2.13(60) 


261 


2.416 




0.5 


0.29 


13.3 


3.64(59) 


197 


2.295 


2.322 






17.8 


7.82(59) 


223 


2.348 




0.1 


0.041 


4.5 


2.06(58) 


127 


2.104 


2.148 






7.3 


7.58(58) 


155 


2.191 




0.05 


0.016 


2.7 


5.32(57) 


103 


2.012 


2.068 






5.0 


2.78(58) 


133 


2.123 




0.01 


0.002 


0.9 


2.69(57) 


64 


1.80 


1.88 






2.1 


2.61(56) 


92 


1.96 





colours refer to the central region of the galaxies (cf. 
Burstein et al. 1988). 

Our model galaxies with AY/AZ=2.5 emit UV radi- 
ation by old H-HB stars only for ages older than about 
5.6 Gyr and in presence of suitable metallicities. Exclud- 
ing the very initial stages with ongoing star formation 
in which strong UV flux can be generated, the colour 
(1550 — V) quickly increases to a maximum and then de- 
creases again (cf. Brcssan et al. 1994, Tantalo et al. 1995). 
This means that in the age range 1 to 6 Gyr the (1550 — V) 
- relation is almost insensitive to the age and very 
close to the 5 Gyr isochrone displayed in Fig. |[ Similar 
considerations hold for the diagram (1550 — V) - [MgFe] 
of Fig. [k]. 

In the (1550 - V) - Up plane (Fig. |), but for two ob- 
jects (M32 and NGC 584), all remaining galaxies have 
data compatible with ages from 8 to about 15 Gyr. The Hp 
index of M32 and NGC 584 is too blue for their (1550 - V) 
colour. The observed (1550 — V) and Ebj of M32, in par- 
ticular, could be matched only going to ages as young as 
about 1 Gyr. However, in such a case the corresponding 
(B — V) colour would be much bluer than the observa- 
tional value of 0.85. 

No such discrepancy is found in the (1550 — V) - 
[MgFe] plane of Fig. [l(] where all galaxies are compatible 
with old ages from 8 to 17 Gyr. M32 in particular seems 
to have an age of about 13 Gyr. 

To improve upon the above analysis we make use of 
three additional observational hints, i.e. the Elg - logS, 
[MgFe] - logS, and (1550 - V) - log £ relations. The data 
refer to the central regions of the galaxies (the Re/8-data 
set of Table |l|) . For the sake of internal homogeneity, we 
adopt the (1550 — V) and logS data of Burstein et al. 
(1988) limited to those objects classified by the authors 
as quiescent (no evidence of active star formation) . 




Fig. 11. The relation between and velocity dispersion log E 
(in Km sec -1 ). The crosses show the uncertainty affecting the 
observational data. Lines of constant age (in Gyr) are shown. 
Finally, the position of M32 is put into evidence 



In order to estimate the velocity dispersion of galactic 
models we made use of the gravitational potential and 
a rough estimate of the radius of the luminous compo- 
nent. The gravitational potential is taken from Bertin et 
al. (1992) and Saglia et al. (1992), 



n = -G(^)(a + ^n' LD ) (14) 
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where G is the gravitational constant and 



1 



* i * q, } Ai 
-t-t 



M32 



Fig. 12. The relation between [MgFe] and velocity disper- 
sion lo g £ (in Km sec -1 ). The meaning of the symbols is as 
in Fig. 



Fig. 13. The relation between (1550 — V) and velocity disper- 
sion log E (in Km sec -1 ). Lines of constant age (in Gyr) are 
shown. The vertical bars indicate the uncertainty in the theo- 
retical log E for the galaxies with different Ml. The cross show 
the uncertainty of the observational data. Finally, the position 
of M32 is shown 



(15) 



For a = 0.5, M L /M D =0.2 and R l /Rd=0.2 as in Bertin et 
al. (1992) and Saglia et al. (1992 ) we get 



Rl 



(16) 



The radius Rl is derived from the Saito (1979a,b) and 
Arimoto & Yoshii (1986, 1987) relationships 



Rt 



26.1 x Mr 



kpc 



(17) 



Finally, the velocity dispersion £ as a function of Ml is 



£ = 260.61 x Ml 



km sec 1 



(18) 



It goes without saying that the theoretical velocity dis- 
persion is affected by a large degree of uncertainty, a rough 
estimate of which can be evaluated considering Ml as a 
sort of parameter ranging from the original value down to 
Ms and taking for £ the mean of the two estimates (£ a d)- 
All these quantities are given in Table |. 

In Fig. [ll] we compare the theoretical - log £ re- 
lations for different values of the age with the data. The 
scatter both in H^ and log £ is large, with galaxies falling 
along isochrones and lines of constant mass spanning an 
ample range of values. However, the hint arises that galax- 
ies with large velocity dispersion are on the average older 
than galaxies with low velocity dispersion. In this diagram, 
M32 has an age of about 6-7 Gyr. 

In Fig. [l^ we compare the theoretical [MgFe] - log £ re- 
lations with the data. But for a general agreement between 
the slopes of the observational and theoretical relation, the 
scatter of the data is so large that firm indications about 
the underlying ages are not possible. For the particular 
case of M32 an old age is most likely. There are galaxies 
(NGC 4489 and NGC 7562) that apparently deviate from 
the mean trend. 

Finally, the theoretical (1550 — V) - log£ relations are 
superposed to the data in Fig. [l3|. Considering all the un- 
certainties, we draw the provisional conclusion that all 
objects of this sample are old. However, looking at the 
particular case of M32, it could be either as old as about 
13 Gyr or as young as about 5 Gyr. 

What we learn from the above analysis is that galaxies 
tend to cluster into two groups. In the first one (NGC 4649 
as prototype) galaxies have old ages (say from 8 to 17 Gyr) 
and normal behaviour in the various planes. We cannot say 
whether the large age range is real or caused by the many 
uncertainties affecting the analysis. In the second group 
(M32 as a prototype), galaxies alternatively change from 
old to young according to the relationship under consid- 
eration, which is an obvious point of contradiction. 
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To explore further the cause of the above discrepancy, 
we look at the three-dimensional Hp - [MgFe] - (1550 — V) 
relation for the small group of galaxies for which the three 
parameters are simultaneously known (11 objects in total, 
cf. Table [I]) , and and compare it with the galactic models 
of 0.5 and 3 M u . The - [MgFe] - (1550 - V) relation is 
shown in Fig. |lj for the Re/8-data. The various lines con- 
nect loci of constant age (thin) and constant mass (thick) . 

The inspection of the data in this space reveals that 
there is a group of galaxies (NGC 4649 as a prototype) 
whose age is confined in the range 13 -r- 15 Gyr (no better 
estimate is possible) and whose colour (1550 — V) is nor- 
mal and compatible with the theoretical expectation for 
old objects containing a certain fraction of high metallicity 
stars. Differences from object to object can be accounted 
for by differences in mass and hence mean metallicity and 
perhaps in age by as much as 3 Gyr (in agreement with 
the estimate from the colour dispersion <7m_v))' 

There is another group of galaxies (M32 and NGC 584 
as prototypes) whose Hp is too blue for their [MgFe] and 
(1550 — V). It seems as if in the Hp - [MgFe] plane an 
old object has been shifted along the line pertinent to 
its mass by some recent episode of star formation which 
has changed without changing [MgFe] and (1550 — V) 
significantly. The episode should not have occurred more 
recently than about 1 Gyr ago, otherwise the (1550 — V) 
colour would have changed. 

Among the various hypotheses that could be invoked 
to explain the puzzling distribution of data in the space of 
the parameters Ha, [MgFe], (1550 — V), and logS, namely 
very low metallicities, large range of absolute ages, and 
burst of star formation, this latter seems to be the only 
viable solution. 

We have already shown in Fig. [?] and Fig. [s] how in the 
Hp - [MgFe] plane an old, red galaxy would be rejuvenated 
by the occurrence of a burst of stellar activity. The exten- 
sion and duration of the loops in the Hp - [MgFe] plane 
depends on the intensity and duration of the burst, and 
whether or not chemical enrichment takes place. Since the 
detailed exploration of space of parameters characterizing 
a burst is beyond the scope of this study (in any case no 
independent observational hints are available), we limit 
ourselves to a few general considerations. 

In the rather idealized case that a burst of star for- 
mation can occur with no accompanying chemical enrich- 
ment, the net effect would be a narrow loop along the 
evolutionary path of a galaxy of given mass. The recovery 
time scales for and colour (1550 — V) are slightly dif- 
ferent: the former depending on the intensity of the burst 
can be as long as 1.5 Gyr, whereas the latter as soon 
as star formation is over quickly goes back to the previous 
values. As chemical enrichment is excluded, no effects (or 
very small) can be seen in the [MgFe] index. 

In the more realistic case, in which chemical enrich- 
ment accompanies the star burst, the evolution in the 
- [MgFe] plane is more complicated, because both rig and 




Fig. 14. The three dimensional correlation between H^, 
[MgFe] and (1550 - V) using the Re/8-data of Gonzales (1993) 
and the UV-fluxes of Burstein et al. (1988). The data are shown 
by the vertical arrows. The thick lines are the loci of constant 
mass (in units of Ml), whereas the thin lines are the loci of 
constant age (in Gyr) of the galactic models. These loci are also 
projected onto the H/j - [MgFe] plane for the sake of clarity. 
See the text for more details 

[MgFe] are now changed by the burst. The increase in 
metallicity would shift the galaxy (or part of it) to redder 
values of [MgFe] while performing the loop in the Hp - 
[MgFe] plane. 

The hypothesis of a burst of star formation while of- 
fering an elegant way out to the difficulties encountered 
above, renders the interpretation of these data more un- 
certain and arbitrary. 

The key difficulty with this hypothesis, is that in the 
case of a primordial origin of elliptical galaxies (all coeval 
within the age range indicated by the CMR), a mechanism 
synchronizing the bursting activity among galaxies would 
be required. 

In the other alternative that all or a number of ellip- 
tical galaxies are the product of a merger, the star burst 
hypothesis is viable and can account for the scatter along 
the Hp axis. The only marginal requirement to be satisfied 
is that the bursting activity does not significantly change 
the global pre-burst metallicity distribution (see the nar- 
row range of metallicities implied by the distribution in 
the Hp - [MgFe] plane). This is the case in which most of 
the stars in the merging galaxies are formed at an early 
epoch and some small amount of star formation occurs 
at the merger epoch. The opposite alternative, in which 
most of the stars are formed in the burst, while easily 
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explaining the scatter in H^, could perhaps introduce a 
scatter in metallicity much larger than that inferred from 
the observations. 

8. Age and metallicity gradients in galaxies 

The galaxies in our sample possess gradients in and 
[MgFe] passing from the Re/2 to the Re/8-data set (cf. 
Gonzales 1993 and the entries of Table [j]) that likely re- 
flect gradients in age and/or metallicity. In order to clarify 
the role played by these gradients in understanding the 
structure and the evolution of galaxies, we perform the 
following analysis. Let us define the quantities 

AlogH / g w = logH^ JV -logH^ w (19) 



A log [MgFe] 



NW 



log [MgFe] N - log [MgFe] w 



(20) 



Fig. 15. The AlogH^^ versus A log [MgFe] NW relation. 
The two arrows centered on (0,0) indicate the ageing and en- 
riching vectors of galaxies and fix the system of coordinates 
used to convert AH^^ and AfMgFe]^,^ into AtNW and 
AZnw, respectively. The cross shows the estimated uncer- 
tainty of the data. See the text for more details 



Fig. 16. The AtNW versus AZnw relation. The meaning of 
the four quadrant's is discussed in the text. The solid line is 
the linear best-fit of the data 



where N and W stand for the Re/8- and R/2-data set, i.e. 
for the central regions and the whole galaxy, respectively. 

The relation between AlogH^^ and 
A log [MgFe] is shown in Fig. [ll| in which the cross 
shows an estimate of the typical uncertainty affecting the 
data. The size of the error bars is assumed to be twice the 
mean value of the errors in the data of Fig. [j]. 

Despite the above uncertainty, with the aid of equa- 
tions (13) through (16) we draw in Fig. [l5| the system of 
coordinates represented by the two vectors AZ — 0.02 
and At = —5 Gyr and centered on the (0,0) point of the 
AlogH^^ - A log [MgFe]^^ plane. The new system of 
coordinates represents the ageing and enriching vectors 
of galaxies or parts of these. Projecting the data of the 
A\ogKp NW - A log [MgFe] NW plane onto the new sys- 
tem of coordinates we get the mean age and metallicity 
differences, At^w and AZ^w, respectively, between the 
central regions and the whole galaxy. These quantities are 
given in Table [7] and their correlation is plotted in Fig. [l6| 

The plane of Fig. |l^ is divided in four quadrant's char- 
acterized by the sign (either positive or negative) of Ai/vw 
and AZ^w- The meaning of the four quadrant's is as fol- 
lows: 

— At]yw < and AZ^w > 0: the nucleus is younger 
and more metal-rich than the external regions of the 
galaxy. This corresponds to a sort of out-inward pro- 
cess of galaxy formation, in which star formation in the 
nucleus continued for significant periods of time, up 10 
Gyr in some cases. Looking at number of galaxies in 
this regions, this case seems to be the most probable 
in nature. 

— Atww > and AZ^w > 0: the nucleus is expected 
to be older and more metal-rich than the external re- 
gions. Only a few galaxies are found in this quadrant 
but with small differences in the age and metallicity. 
These galaxies could have been formed by a sort of 
in-outward mechanism on a rather short time scale. 
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Table 7. Data on the AlogH^^ - A log [MgFeJ^^, plane and mean metallicity and age differences between the nucleus and 
the whole galaxy, AZnw and Atisrw, respectively. 



NGC 


A log [MgFej ww 


A log Hg ww 




Ai jvw 


NGC 


A log [MgFeJ ww 


A log Hg ww 


AZjvw 






221 


+0.00 


+0.03 


+0.006 


-3.962 


4472 


+0.01 


-0.01 


+0.002 


+0.679 


224 


+0.01 


-0.01 


+0.002 


+0.679 


4478 


+0.03 


+0.02 


+0.015 


-4.566 


315 


+0.03 


-0.01 


+0.009 


-0.604 


4489 


+0.05 


+0.02 


+0.023 


-5.849 


507 


+0.02 


-0.07 


-0.006 


+7.962 


4552 


+0.03 


-0.01 


+0.009 


-0.604 


547 


+0.03 


+0.05 


+0.021 


-8.528 


4649 


+0.03 


+0.01 


+0.013 


-3.245 


584 


+0.03 


+0.01 


+0.013 


-3.245 


4697 


+0.06 


+0.02 


+0.026 


-6.491 


636 


+0.03 


+0.01 


+0.013 


-3.245 


5638 


+0.04 


-0.01 


+0.013 


-1.245 


720 


+0.01 


-0.11 


-0.017 


13.887 


5812 


+0.02 


+0.00 


+0.008 


-1.283 


821 


+0.03 


-0.04 


+0.004 


+3.358 


5813 


+0.02 


+0.06 


+0.019 


-9.208 


1453 


+0.03 


-0.03 


+0.006 


+2.038 


5831 


-0.01 


-0.01 


-0.006 


+ 1.962 


1600 


+0.00 


-0.05 


-0.009 


+6.604 


5846 


+0.04 


+0.06 


+0.026 


-10.491 


1700 


+0.03 


+0.00 


+0.011 


-1.925 


6127 


+0.03 


+0.00 


+0.011 


-1.925 


2300 


+0.02 


+0.02 


+0.011 


-3.925 


6702 


+0.02 


-0.01 


+0.006 


+0.038 


2778 


+0.03 


+0.06 


+0.023 


-9.849 


6703 


+0.04 


+0.01 


+0.017 


-3.887 


3377 


+0.06 


-0.01 


+0.021 


-2.528 


7052 


+0.03 


-0.08 


-0.004 


+8.642 


3379 


+0.03 


+0.01 


+0.013 


-3.245 


7454 


+0.05 


+0.01 


+0.021 


-4.528 


3608 


+0.04 


-0.01 


+0.013 


-1.245 


7562 


-0.03 


-0.01 


-0.013 


+3.245 


3818 


+0.06 


-0.03 


+0.017 


+0.113 


7619 


+0.03 


-0.04 


+0.004 


+3.358 


4261 


+0.03 


+0.02 


+0.015 


-4.566 


7626 


+0.04 


+0.00 


+0.015 


-2.566 


4278 


+0.03 


-0.03 


+0.006 


+2.038 


7785 


+0.02 


+0.03 


+0.013 


-5.245 


4374 


+0.02 


-0.01 


+0.006 


+0.038 













— AtNw > and AZnw < 0: the nucleus is older but 
less metal- rich than the external regions. These ob- 
jects do not find a straightforward explanation. They 
could correspond to cases of merge in which an old 
compact primeval object with scarce chemical enrich- 
ment subsequently captured other smaller galaxies in 
which chemical enrichment has already proceeded to 
substantial levels. 

— AtNw < and AZnw < 0: the nucleus should be 
younger and less metal-rich than the external regions. 
No galaxy is found in this quadrant. This kind of 
galaxy structure and star formation in turn is not al- 
lowed. 

— There appears to be a tight correlation between 
AZnw and AtNw as indicated by the linear best-fit 

AZnw = -0.0018 x At NW + 0.0073 (21) 

with AtNw m Gyr. This relation indicates that galax- 
ies as a whole or their constituent parts in spite of 
the large variety of properties (masses, dimensions, 
colours, luminosities, etc.) are ultimately governed by 
age and metallicity, and that chemical enrichment in 
these obeys an universal law. 

This analysis cannot yield the true ages and metallic- 
ities of the nuclear and peripheral regions of a galaxy but 
only their mean difference. 

It goes without saying that the above results de- 
pend on the transformation from the AlogH/3 NW - 
A log [MgFe]^^ to the AtNw ~ AZ N w plane. However, 
the transformation in use cannot be grossly in error be- 
cause it stems from the properties of SSPs and model 
galaxies. As already amply discussed in the previous sec- 
tions both match the data in the basic Ha - [MgFe]. 



The numerous group of galaxies in the first quadrant 
(AtNw < and AZnw > 0) is of particular interest. In 
spite of the afore mentioned uncertainty, for these galaxies 
we suggest that a sort of out-inward process of formation 
ought to have occurred, lasting in many cases up to sev- 
eral Gyr. In this context, it would be interesting to know 
whether the time AtNw correlates with other important 
parameters such as the velocity, dispersion log E, the mass 
to blue luminosity ratio (M/Lb)q, the effective radius Re, 
and the total mass. The inspection of the data contained 
in Tables @ and @, indicates that there could be a weak 
correlation between AtNw and the velocity dispersion E: 

At NW = 3.327 x log E - 10.65 (22) 

where AtNw is in Gyr, E is in km sec -1 , and the corre- 
lation coefficient is 0.43. This relation has been derived 
discarding the four objects with the longest, perhaps sus- 
pected, At NW , namely NGC 547, NGC 2778, NGC 5813 
and NGC 5846. Similar but weaker dependencies on the 
total mass and effective radius can be found. 

Therefore, the hint arises that the time difference be- 
tween the global and the nuclear activity decreases at in- 
creasing velocity dispersion and perhaps total mass. This 
time difference can be taken as an indication of the overall 
duration of the star forming period. 

9. Concluding remarks 

In this paper we have addressed the question whether age 
and metallicity effects on the stellar content of elliptical 
galaxies can be disentangled by means of broad-band col- 
ors and line strength indices derived from detailed chemo- 
spectro-photomctric models. The main goal of our study 
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is to cast light on a number of interwoven questions: are 
elliptical galaxies made only of old stars generated in a 
dominant initial episode of star formation? is there any ev- 
idence of subsequent star forming episodes ? If so, has this 
activity an internal origin or is it induced by the merger 
mechanism ? 

The analysis of the observational data in the space of 
the parameters Hp, [MgFe], logS, and (1550 — V) and of 
the gradients inH^ and [MgFe] within individual galaxies 
allow us to draw a number of conclusions. 

— The distribution of galaxies in the Hp - [MgFe] plane 
does not correspond to a mere sequence of metallic- 
ity with bluer galaxies significantly more metal-poor 
than the red ones. Metallicity effects are expected to 
broaden the distribution mainly along the [MgFe] axis. 

— Equally, the distribution of galaxies in the Hp - [MgFe] 
plane does not constitute a mere sequence of age with 
bluer galaxies significantly younger than the red ones. 
Although some scatter in the age is possible, it is not 
likely to be as large as indicated by the formal match of 
data with isochrones in this diagram. M32 in particu- 
lar according to its Hp should have an age of about 
3 Gyr, in disagreement with its (1550 — V)=4.5 for 
which much older ages are needed according to the 
present theory of population synthesis. 

— The observed distribution of galaxies in the H^ - 
[MgFe] plane does not agree with the CMR-strip on 
the notion that they evolve as passive objects in which 
star formation took place in an initial episode. Most 
likely, the history of star formation was more compli- 
cated than this simple scheme. 

— The observed (1550 — V) colours are not compatible 
with ages younger that about 6 Gyr. UV excesses 
stronger than (1550 — V)=5.5 at younger ages are pos- 
sible only if the age is younger that 1 — 2 Gyr. However, 
in such a case the remaining UBV colours would be 
too blue compared with the observational data. Most 
likely, all galaxies in the sample are globally older than 
at least 6 Gyr. 

— As far as ages are concerned, our analysis is still un- 
able to give a definitive answer. Besides the group 
of galaxies like M32 and NGC 584 whose properties 
lead to contrasting results depending on the particular 
relationship among the four parameters Hg, [MgFe], 
(1550 — V) and S under examination, the remaining 
galaxies appear to span a wide range of ages going 
from 8 to 15 Gyr. However, this large scatter can be 
partly due to observational uncertainties and partly 
to a real age spread. It is worth pointing out that for 
the very small subgroup of galaxies with all the four 
parameters simultaneously available, the objects with 
an old age (13-^15 Gyr) like NGC 4649 show an age 
spread amounting only to about 3 Gyr, which is con- 
sistent with the small scatter in the CMR of cluster 
galaxies. 



— The present analysis cannot yet cast light on the 
dominant mechanism by which elliptical galaxies are 
formed. However, the following considerations can be 
made. If elliptical galaxies are all old, coeval and evolv- 
ing in isolation, the difficulties encountered with the in- 
terpretation of the Ha, [MgFe], (1550 - V), and logS 
data can be removed by invoking later mild episodes of 
star formation. The major problem with this sugges- 
tion is that in the case of evolution in isolation, some 
mechanism synchronizing the bursts of star formation 
should be required in order to dislocate galaxies from 
the CMR-strip to their observed location in the H^ - 
[MgFe] plane. Indeed, even if the statistics is limited, 
no galaxies are seen along the CMR-strip but for those 
at the red end. A simple way out could be that the 
star forming activity is limited to the nuclear region, 
thus affecting the Hp index without changing signifi- 
cantly the global value of the broad-band colours (cf. 
the case of NGC 4374 and NGC 4697 for which little 
age scatter is indicated by the CMR, whereas a signif- 
icant age scatter is suggested by the indices). No such 
difficulty would exist if already formed ellipticals suf- 
fer in more recent times from interaction with other 
small size galaxies (for instance the capture of small 
gas rich systems) thus inducing some mild star forma- 
tion activity with some chemical enrichment. Equally, 
if elliptical galaxies are the result of classical mergers of 
spiral galaxies. In such spread in Hp mimicking 
an age sequence for the resulting composite popula- 
tion would naturally follow. The accompanying chem- 
ical enrichment can be easily adjusted to match the 
small observational scatter along the [MgFe] axis. 

— Although based on a handful of objects, there is a in- 
teresting scenario emerging from the analysis of the 
four dimensional space H^, [MgFe], (1550 — V), and 
log X. In brief, galaxies of high velocity dispersion tend 
to be confined at the red end of the Hg - [MgFe] distri- 
bution. Their properties are consistent with being old 
objects (some scatter in the age is perhaps possible) 
that are able to exhaust the star forming process at 
very early epochs with no need for later episodes of 
stellar activity. In contrast, galaxies of lower velocity 
dispersion have more contrasting properties. It seems 
as if in these systems star formation, following the ini- 
tial activity most likely at epochs as old as in other 
galaxies, later went through a series of episodes tak- 
ing place in different epochs that vary from galaxy to 
galaxy. Equivalently, it could be that star formation 
continued perhaps at minimal levels over significantly 
longer periods of time. This would mimic a sort of age 
sequence. 

— The analysis of the gradients in Hp and [MgFe] within 
individual galaxies has revealed that for the majority 
of these the nucleus is younger and more metal-rich 
than the peripheral regions, thus resembling a out- 
inward mechanism of galaxy formation, and that the 
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global process of star formation may last very long, 
with duration varying from galaxy to galaxy. It seems 
that the duration of the star formation activity is in- 
versely proportional to the galactic velocity dispersion 
and perhaps mass. There are two less numerous sub- 
groups for which the nucleus is older and more metal- 
rich or older and less-metal rich than the periphery of 
the galaxy. The interpretation of these structure is not 
straightforward and it is left to future investigation. 

— The fact that most galaxies are found in the group 
characterized by a younger and more metal-rich nu- 
cleus calls to mind the fact that recent bursts of star 
formation in the merger scenario are expected to peak 
in the centers of elliptical galaxies (Alfcnslcbcn & Ger- 
hard 1994). However, an internal cause related to the 
building up process of galaxy formation cannot be ex- 
cluded as perhaps suggest by the correlation between 
velocity dispersion and ISi^w- 

— The idea that the overall duration of the star forming 
activity is inversely proportional to the velocity disper- 
sion and perhaps mass of galaxies does not contradict 
the current information of abundance ratios in ellipti- 
cal galaxies (Carollo et al. 1993, Carollo & Danziger 
1994, Matteucci 1994). Indeed it would be consis- 
tent with the expected trend for the ratio [Mg/Fe] 
inferred from the narrow band indices (Faber et al. 
1992; Worthey et al. 1992; Davies et al. 1993). From 
the comparison of the observed indices with model 
indices, assuming solar partition of elemental abun- 
dances, it is concluded that the average [Mg/Fe] in 
giant elliptical galaxies exceeds that of the most metal 
rich stars in the solar vicinity by about 0.2-0.3 dex. 
Furthermore, the ratio [Mg/Fe] is expected to increase 
with the galactic mass up to the this value. Accord- 
ing to Matteucci (1994) a possible explanation of this 
trend is that the efficiency of star formation is an in- 
creasing function of galactic mass. As a consequence of 
this, massive ellipticals should have formed on shorter 
time scales than smaller ellipticals. In brief, the con- 
stant super-solar value for giant ellipticals hints that 
only an unique source of nucleosynthesis is contribut- 
ing to chemical enrichment. This is just the case of 
massive stars exploding as type II supernovae and thus 
releasing elemental species (Mg and Fe in particular) 
in fixed proportions. This implies a short time scale 
of star formation, i.e. of about 0.1 Gyr. The trend de- 
creasing with galactic mass means that another source 
of Fe intervenes altering the [Mg/Fe] ratio. It is worth 
recalling that Mg is produced by massive stars. There- 
fore, it is the contribution to Fe that must change. 
The goal is achieved invoking type I supernovae gener- 
ated by mass accreting white dwarfs in binary systems. 
The minimum time scale for their formation is longer 
than 0.1 Gyr. This implies a longer duration of the star 
forming period. In this context, we would like to recall 
that in order to match the CMR with our infall mod- 



els we had to increase the efficiency parameter v by a 
factor of 12 as the galactic mass increases from 0.01 to 
3 M12. Furthermore, the duration of the star forming 
activity determined by the onset of galactic winds was 
found to decrease with the galactic mass (see t gw in 
Table ||) . The trend is small but in the right direction. 
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Table 4. Evolution of the integrated indices for SSP. Ages are in Gyr. 
Age Hp [MgFe] < Fe > Mgb Mg2 Mgl Fe 527Q Fe. 
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Table 4. continued 
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